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Abstract 

We present the first numerical implementation of a non-perturbative renormal- 
ization method for lattice operators, based on the study of correlation functions in 
coordinate space at short Euclidean distance. The method is applied to compute 
the renormalization constants of bilinear quark operators for the non-perturbative 
0(a)-improved Wilson action in the quenched approximation. The matching with 
perturbative schemes, such as MS, is computed at the next-to-leading order in con- 
tinuum perturbation theory. A feasibility study of this technique with Neuberger 
fermions is also presented. 



1 Introduction 

Correlation functions of quantum-field operators are computed non-perturbatively on the 
lattice by Monte Carlo simulations. From their long-distance behaviour, in QCD, non- 
perturbative features of the underlying theory, such as the hadron spectrum and matrix 
elements, can be extracted. On the other hand, their behaviour at short distance is ex- 
pected to be controlled by perturbation theory and by the operator product expansion 
(OPE). As a result, relevant ultraviolet informations of the theory, such as the range of 
applicability of the OPE, the values of the Wilson coefficients and the renormalization 
constants of composite operators, can in principle be obtained by a comparison of the 
perturbative formulae with the numerical results at short distance. This is the basic idea 
behind the non-perturbative renormalization techniques proposed in the last decade and 
widely used in present simulations. With the RI/MOM method, renormalization condi- 
tions are imposed on quark and gluon Green functions computed non-perturbatively in 
momentum space, in a fixed gauge, with given off-shell external states of large virtual- 
ity pP-jni- In the Schrodinger functional (SF) approach, the renormalization conditions 
are imposed in coordinate space at a given finite physical distance on suitable gauge in- 
variant correlation functions in finite volume with SF boundary conditions JI]-^!]. The 
step-scaling technique ^2| can then be used to convert the renormalization constants to 
their renormalization group invariant definitions. An advantage of these methods is that 
the conversion to more popular continuum schemes, such as the MS scheme, can be imple- 
mented by performing a calculation only in continuum perturbation theory, by comparing 
renormalized correlation functions at short distances computed in dimensional regulariza- 
tion in the two schemes. In this way, more tedious calculations with lattice perturbation 
theory are completely avoided. 

In this paper we present the first numerical implementation of a non-perturbative renor- 
malization method based on the study of lattice correlation functions at short Euclidean 
distances in coordinate space ^Tj. We call this approach the "X-s]9ace" scheme. Prelimi- 
nary results were presented in Refs. fHlCSl- We have computed numerically the two-point 
functions of all dimension three bilinear quark operators by discretizing the gluons a Id Wil- 
son and the fermions with the non-perturbatively 0(a)-improved Wilson action. A feasibil- 
ity study of this technique with Neuberger fermions is also presented. The multiplicatively 
renormalization constants of lattice bilinear operators are evaluated non-perturbatively by 
imposing renormalization conditions directly in x-space at short distance. The condition 
Xq ^ a, where a is the lattice spacing and Xq is the renormalization point, has to be satis- 
fied in order to keep discretization errors under control. On the other hand, the matching 
of the renormalization constants to the MS scheme (or any other continuum scheme) can 
be computed in continuum perturbation theory when xq <C Aq^j-,. In this study we show 
that for a < 0.05 fm (i.e. 1/a > 4 GeV) it is possible to find a region on the same lat- 
tice in which perturbation theory can be applied and discretization effects are still under 
control. The existence of the window a < Xq < AqJ^j-, requires however rather fine lattices, 
large volumes and therefore expensive simulations. Alternatively, one can appeal to a step- 
scaling technique analogous to the one proposed in Ref. ^21- The X-space renormalization 
method involves only gauge-invariant correlation functions among local operators at finite 



physical distance, and can be easily applied to any fermion discretization. It can be very 
powerful for the evaluation of the renormalization constants of composite operators, such 
as the four-fermion operators relevant for the phenomenology of hadronic weak decays. 
The matching to the MS scheme can be performed by using only continuum perturbation 
theory and the method is very simple to implement. 

The paper is organized as follows: in sec. |21 we summarize the relevant formulae in 
continuum perturbation theory and define the renormalization conditions which we will 
use in coordinate space; in sec. El we discuss the numerical results and in sec. El we present 
our conclusions. 

2 X-space renormalization and perturbation theory 

In this section we define the X-space renormalization scheme for bilinear quark opera- 
tors and provide the perturbative expressions, at the next-to-leading order (NLO), needed 
to convert the results to any other continuum renormalization scheme, such as the MS 
scheme^. 

We consider the correlation functions of flavor non-singlet bilinear quark operators of 
the form 

(Or(x)Or(O)) , (1) 

where 

Orix) = ^(x)r^(x) (2) 

with Or = {S,P,V^,Af,Tf,^} for T = {1, 75, 7^, 7^75, |[7^,7i^]} respectively and with 
flavor indices omitted. 

Following P2j, we impose non-perturbatively, in x-space and in the chiral limit, the 
renormalization conditions 



lim(O-(x)O-(0)) ^ = (Or(a;o)Or(0))rt (3) 

where the renormalized operator is Op (a;, xq) = Z^ixo) Or(x) and xq is the renormalization 
point. The renormalization condition ^ defines the X-space scheme. Note that xq must 
satisfy the condition a <^ xq <^ ^qcd ^'^ keep non-perturbative and discretization effects 
under control. 

In order to illustrate the procedure and get the expressions needed to convert to the 
more popular continuum schemes, we have computed the correlation functions, at two loop 
in naive dimensional regularization (NDR), in the massless case. The results, in Euclidean 

^AU formulas presented in this section are obtained in the infinite- volume limit. When correlation 
functions in small volumes are considered, their perturbative expressions may need to be modified according 
to the boundary conditions used ^J E] 



space, read 
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{-I^HuyXjlpay^)) ^ 2\3 1 2 Mi^PO" ^ ^lupa j 



where 



1 + + 4 - - 7i5 - ^^ ln(^- — 

47r I 3e 3 ' 2 W 



■^ 111/ pa ~ Op^pOi/a dpaOup 



rp(2) ■^p-^pr ■^P'^fT r •'^U-'^pr ■^U'^a r /p.\ 

-'■ uupa ~ ~~3~ '^" Z^ '^P Z^ P-" Tj^ P-P ' ^ ) 



IM/pa 2 "'' ^2 "P ^2 P'^ • ^2 

tAj lAj tXj iXj 



In these expressions a, is the strong couphng constant, 1/e = 1/e — ln(47r) — 7^; (we define 
d = A — 2e the space-time dimension) and the LO anomalous dimensions are 7^ = —8, 
7ji = 8/3, and 7^ = 0, the latter due to the conservation of the vector current. This 
conservation also determines the tensor structure of the vector current correlator. Note 
that, as expected, the leading short distance behaviour of the correlation functions in 
Eq. (JH) is governed by (x^)~^. For the scalar and vector correlators the results in Eqs. (J3J) 
agree with previous computations [201 HH 122] ■ 

By imposing the renormalization conditions to the results in Eq. (@)), we obtain the 
correlation functions for the renormalized operators in the X-space scheme: 

(5^(x,Xo)5^(0,Xo)) = -^^K^ix^x,) 

2Nc (I XpXu 



(V;^(x,Xo)K^(0,Xo)) = --±^^1^^^^^ (6) 

(r^(x,xo),.r^(o,xo)p.) = -^^^Qr«.-rS.)i^^(x,. 



with 



Af(...„)^l -of £..>(!). (7) 



In the MS scheme, on the other hand, the correlation functions of the renormalized 
composite operators, 0^^(a;, /i) = Z^^(/i) Or(x), are obtained by subtracting the pole 1/e 
on the r.h.s of Eq. Q. The relations between the MS and the X-space scheme are thus 
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This also implies that the renormalization constants in the two schemes are related by 
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The renormalization condition Q does not satisfy the vector and axial vector Ward Iden- 
tities. At the NLO this can be easily seen by noticing that in the MS scheme, which 
preserves them, Ky'^{xo,fi) has a finite term proportional to a,. In the X-space scheme 
this contribution is included in the renormalization constant and therefore Zy^^/Zy 7^ 1, 
i.e. the Ward Identities are broken. They can be recovered by using continuum per- 
turbation theory, which for the vector correlator is known up to four loops j^U 122] 5 or 
non-perturbatively by matching the result for Zy in the X-space scheme with the Ward 
Identity determination |23j . 

We conclude this section by recalling the expression for the renormalization group 
evolution of the renormalization constants at the NLO in a,: 



z^'in') 



„MS 



(/^O 



.MS 



(/^) 



zn^^) 



where 



„(0) 



MS 



(/i)=a,(/i)2/'o 



Att 




(0)' 






(10) 



(11) 



with as defined in MS scheme. In the quenched theory, i.e. with Nf - 
two coefficients of the expansion of the /3-function are /3o = 11 and Pi 
the two-loop anomalous dimensions in the MS scheme are 7^ = —404/3, 



0, the first 
102, while 

L' 
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'■^^ - and 



7j, = 724/9 j21l 123 ESI • In our study, the matching between X-space and MS schemes has 
been performed at a scale /i ~ I/xq. With this choice, logs in Eqs. © are small and need 
not to be resummed. In contrast, larger logs enter the evolution from the scale ^ ~ I/xq 
to the conventional scale /z = 2 GeV at which our final results are quoted. For this reason, 
the evolution function in Eq. (fTTj) has been resummed by using the renormalization group 
equations at the NLO. 
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Figure 1: The vector correlator Cvv{x) for different values ofx^ as a function oftlie quark mass. 
The line represents the result of the linear extrapolation to the chiral limit. 

3 Numerical results 

In this section we provide the numerical details of our computation and present the re- 
sults obtained non-perturbatively for the renormalization constants of the bilinear quark 
operators. 

We generated a sample of 180 gauge configurations in the quenched approximation with 
the standard SU(3) Wilson gluonic action at /5 = 6.45 (a ~ 0.048 fm) and y = 32^ x 70. For 
these configurations we evaluated fermion propagators with the non-perturbatively 0(a)- 
improved Wilson action for hopping parameter values /€=0.1349, 0.1351, 0.1352, 0.1353. 
We computed the two-point flavor non-singlet correlation functions 



Cssix) = {S{x)Sm , 
Cvvix) = J2(^,{x)V,m 



Cpp{x) = {P{x)P{0)) , 

Caa{x) = Y,{M^)Mo)) 



(12) 
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of local bilinear operators in the standard way and estimated the statistical errors by a 
jackknife procedure. These functions have been averaged over points which are equivalent 
under hypercubic rotations. In the range of x'^ we have studied, our data show a mild mass 
dependence, and a linear or quadratic extrapolations to the chiral limit give compatible 
results within the errors. In the following we will show the results linearly extrapolated to 
the chiral limit. An example of this extrapolation, in the case of the vector correlator, is 
illustrated in Fig. [H 

In Fig. 121 (left) we show the correlation function Cvv{x) as a function of x^ and the 
corresponding one parameter fit to A {x"^) . The naive expected behaviour is clearly 
satisfied, even if a large scattering of the points due to lattice artifacts is visible, particularly 
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Figure 2: Cvv{x) in the interacting (left) and in the free (right) theory. In the first case, the 
dotted curve is a one parameter fit showing the (x^) behaviour. In the free case the curve 
represents the prediction of the free theory in the continuum hmit 
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Figure 3: The corrected vector correlator Cyy{x). The dotted curve is a one parameter fit 
showing the (x^) behaviour. 



in the short distance region. Lorentz invariance requires the correlator in the continuum 
hmit to be a function of x^ only. The lattice data presented in Fig. |21 (left) show instead 
that the results for Cvv{x) computed at points which correspond to the same value of x"^ 
are often quite separated. To clarify the origin of these effects, we studied the correlators 
in the free theory as a function of volume, lattice spacing and quark masses. The free 
theory prediction for the lattice correlation function Cwi^), in infinite volume and in the 
chiral limit, is shown as an example in Fig. |21 (right). The spread of the data observed in 
the interacting case turns out to be well reproduced in the free theory, at fixed volume and 
lattice spacing. For values of x^ in the perturbative region, finite volume effects are found 
to be negligible in the range of masses we use. On the other hand, the spread of the data 
in the free case is considerable reduced by decreasing the lattice spacing. This suggests 
that the dominant contributions due to lattice artifacts comes from discretization effects. 
In order to reduce discretization effects in the interacting case, we define "corrected" 
correlation functions 

^-(^) - Wj <^^* 

where Arr(a;) is the ratio of free correlator on the lattice over the continuum one, computed 
in infinite volume and in the chiral limit, 



(Or(x)Or(O)) 



free 
lat 



^Mx)= ' ''I 1 - (14) 
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free 
cont 



By construction, Apr (a;) is equal to unit up to discretization effects. The results for the 
corrected function Cyyl^x) are shown in Fig. El These results, as well as those used in 
the following analysis, have been also averaged over points which correspond to the same 
x^. The comparison between Cwi^) and Cyy^x) shows that, once tree-level discretization 
effects are removed, the spread of the data is greatly reduced. This analysis further supports 
the interpretation that the spread in the interacting theory is due to discretization effects. 
Similar conclusions apply also to other correlators. 

Once the correlators have been corrected at tree level with the factor that attempt to 
reduce discretization effects at the leading order, they still suffer for remaining discretiza- 
tion errors, at Olg"^ a^) and higher. Their effects is shown in Fig. HJ where we plot the 
corrected correlation functions Cfp(x) normalized to their continuum counterparts in the 
free theory, 

«rr(x) = -^j!^ (15) 

in the case of the vector current and the scalar density operators. Note that, particularly 
in the scalar case, the remaining lattice artifacts are still larger than statistical errors. A 
further reduction of these effects could be obtained by computing the 0{g^ a?) terms in 
lattice perturbation theory. 

In order to extract the renormalization constants, the previous analysis suggests to 
implement the renormalization condition defined in Eq. Q directly to the corrected cor- 
relation functions Cfp(x). This is equivalent to impose 

(O?(x)O^(0)) 
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Figure 4: The ratios Rvvix) (left) and Rss^x) (right) dehned in Eq. if73)). 

Zy Za Zx 

This work 0.801(2)(18)(6) 0.833(2)(27)(6) 0.895(2)(21) 

RI-MOM 0.803(3) 0.833(3) 0.898(6) 

SF 0.808(1) 0.825(8) 



Zp/Zs 



This work 0.888(2) (8) 
RI-MOM 0.897(4) 
SF 0.912(9) 



0.702(4)(27)(23) 0.624(3)(19)(21) 
0.679(8) 0.609(8) 

0.61(1) 



Table 1: Wilson results for the renormalization constants in the MS scheme at /x = 2 GeV, and 
comparison with other non-perturbative techniques: RI-MOM j^ and Schrodinger Functional 

I23l23ESf- 

1X 



Since 0^{x) = Z^{xo)Or{x), the above condition imphes 



Z^ixo) 



l/^Rvv{xo) 



(17) 



The renormahzation constants Zy-, Za-, Zs, Zp, Zt in the X-space scheme and the 
ratio Zp/Zs are shown in Fig. El as a function of (xo/a)'^. The scattering is reduced 
for [xo/a)"^ > 8, signahng that discretization effects are moderate above this point. For 
(xo/a)^ < 21, which corresponds to I/xq > 0.9 GeV, the dependence of the renormal- 
ization constants on the renormalization scale is compatible with the NLO prediction of 
perturbation theory, indicated by solid lines in Fig. El In particular, the ratio Zp/Zs 
is expected to be independent of the renormalization scale and we see that in the range 
(xo/aY = [9, 21] a reasonable good plateau is observed. In this window we extract the val- 
ues of all renormalization constants by fitting the corresponding correlation functions with 
the perturbative formulae given in sec. |21 The results, translated to the MS scheme at the 
scale yU = 2 GeV are presented in Tab.JH The first quoted error is statistical and it is by far 
the smallest one. The second is an estimate of the uncertainty coming from the spread of 
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Figure 5: Renormalization constants in the X-space scheme as a function of the renorniahzation 
scale xq. The sohd line indicate the scale dependence predicted by NLO perturbation theory. 
The dashed lines show the boundaries of the fitting window . 

the points within the fitting window. This error could be further reduced by going to finer 
lattices or by evaluating the 0{g'^ a^) terms in lattice perturbation theory. The third error 
in Tab. n is an estimate of the systematics due to higher orders in continuum perturbation 
theory, obtained by varying the renormalization scale /i in the perturbative expressions of 
Eq. © in the range 1 < yUXo < 2. This uncertainty can be reduced by performing a N^LO 
computation in perturbation theory and/or by implementing the step scaling technique^ 
proposed in Ref. [^. The latter would require simulations at several lattice spacings and 
goes beyond the scope of this exploratory study. 

In order to investigate the applicability of the X-space renormalization method to dif- 
ferent discretizations of the fermionic action, we also performed a feasibility study by using 
Neuberger fermions. We used 80 configurations generated with the same gluonic action 
at /3 = 6.0 (a ~ 0.093 fm) and V = 16^ x 32 which were retrieved from the repository 
at the "Gauge Connection" ^7\. For these configurations, overlap propagators at bare 
masses ma = 0.040, 0.055, 0.070, 0.085, 0.100 have been evaluated, see Refs. 13 E] for 
details. The results have been quadratically extrapolated to the chiral limit, as suggested 
by perturbation theory. The results for the renormalization constants of the vector and 
axial-vector currents are shown in Fig. IHlas a function of {xo/ay. Since in this case the 
lattice spacing is larger than in the case of Wilson fermions, the window contains at most 
three points, which makes it difficult a reliable comparison with perturbation theory. Nev- 
ertheless the data plotted in Fig. IB] show a smooth behaviour compatible with the chiral 
properties of Neuberger fermions, and the value of Zy = Za for 3 < (xo/a)^ < 5, corrected 
for the matching factor in Eq. Q, is compatible with the Ward identity estimate obtained 

^It is interesting to note that, since the running of the operators with the renormaHzation scale is 
scheme dependent but rcgularization independent, it is possible to implement the step scaling technique 
by using any discretization of the fermionic action. The results, extrapolated to the continuum limit, can 
then be used to evolve the renormalization constants computed non-perturbatively also with a different 
discretization. 




Figure 6: Renormalization constants Zy, Za and the ratio Zy JZa for Neuberger fermions as a 
function of xq. 

in Ref. [3 El, Za = 1.55(4). In this case, the implementation of a step scahng technique 
would have allowed us to impose the renormalization conditions at shorter distances, where 
perturbation theory is reliable but discretization effects remain negligible. 

4 Conclusions 

We have studied the correlation functions of two fermion bilinear operators in coordinate 
space at short Euclidean distances. A good statistical signal can be obtained with a small 
number of configurations. The spread of the data computed at different lattice points with 
the same x"^, indicates that discretization errors can be large. This spread has been greatly 
reduced by normalizing the correlation functions with the analogous ones computed in the 
free theory at finite lattice spacing. A straightforward application of these results is a 
determination of the renormalization constants of the composite operators. Even if with 
larger uncertainties, the values of the renormalization constants which we have obtained 
are in good agreement with previous non-perturbative determinations. This technique can 
be easily applied to any fermion discretization, as shown in this paper, and it involves only 
gauge-invariant correlation functions among local operators at finite Euclidean distance. 
The matching with more popular renormalization schemes, such as the MS scheme, is 
easy because it requires calculations only performed in continuum perturbation theory. 
Therefore in the future, after more accurate studies of the systematic, the X-space method 
could become a powerful technique to renormalize the four-fermion operators relevant in 
hadronic weak decays. 
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